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SUMMARY 


This  report  is  concerned  with  special  problems  which  arise  in 
connection  with  calculations  of  the  transport  properties  of  partially 
ionized  monatomic  gases.  The  most  serious  problem  concerns  the  lack  of 
agreement  of  the  usual  thermal  conductivity  expression  in  the  limit  of 
full-ionization  with  other  results  derived  explicitly  for  this  case. 

It  is  shown  that  satisfactory  agreement  can  be  obtained  in  this  limit  if 
one  uses  the  third  rather  than  the  second  approximation  in  the  Chapman- 
Enskog  theory .  Expressions  for  the  fourth  and  lo*  *  approximations  to 
the  thermal  conductivity,  the  thermal  diffusion  c  fficier.t,  and  the 
ordinary  diffusion  coefficient  of  multicomponent  gases  are  derived. 

The  viscosity  of  this  mixture  is  considered  to  the  second  approximation. 
It  is  shown  that  the  thermal  diffusion  plays  a  veiy  importar-.":  role  in  an 
ionized  gas.  Neglect  of  this  effect  can  cause  the  therma’  crnductivity 
to  be  seriously  overestimated.  An  expression  is  presented  %hich  approxi¬ 
mates  the  effect  of  the  thermal  diffusion  on  the  thermal  cc r.ductivity. 

The  charged  particle  cross-section  is  considered  with  the  screened 
Coulomb  potential.  Cross-sections  derived  with  this  potential  should  be 
acc\irate  for  denser  plasmas  that  can  be  accurately  considered  with  the 
usual  Debye  cut-off  cross-section.  Convergence  of  the  higher  approxi¬ 
mation  is  checked  by  applying  the  formulai.  to  the  calculation  of  the 
properties  of  several  special  mixtures  where  the  molecules  obey  the 
Coulomb  and  other  inverse-power  interaction  potentials. 
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1.  INTRODUCTION 


In  the  past  sevei’al  years  there  has  been  increased  reliance  on 
computed  rather  than  measured  values  of  transport  properties  of  gases. 

This  has  come  about  because  of  difficulties  in  measuring  these  coeffi¬ 
cients  with  any  accuiacy  at  the  higher  temperatures  presently  of  interest. 
The  ceLlculations  start  from  potentials  of  interaction  between  molecules, 
usually  determined  from  molecular  beam  scattering  experiments,  and  use 
theoretical  or  approximate  formulas  for  the  computations.  The  theoreti¬ 
cal  methods  are  usually  the  first  or  second  approximations  in  the  Chapman- 
Enskog  theory  [Chapman  and  Cowling,  1952;  Hirschfelder,  Curtiss  and  Bird, 
I96I1]  .  The  approximate  formulas  [Wilke,  1950:  Mason  and  Saxena,  1958: 
Brokaw,  1958]  are  usually  derived  from  the  theoretical  formulas  by  making 
certain  approxl.nations.  The  approximations  are  ultimately  justified  by 
comparison  of  calculations  made  with  the  resulting  formulas  with  those 
made  with  the  theoretical  formulas.  The  agreement  is  fairly  good  for 
most  low  temperat’ore  gases. 

Both  the  theoretical  and  approxLmate  formuL.s  have  been  applied 
recently  to  the  calculation  of  properties  of  dissocia'cing  hydrogen 
[Vandersllce  et  al,  I962]  and  nitrogen  and  oxygen  [Yun  et  al,  I962] . 

Later  these  formulas  have  been  applied  to  gases  in  which  ionisation  is 
also  important.  Eastlund  [1983]  calculated  the  transport  procerties  of 
partly  ionized  cesium  vapor  at  temperatures  up  to  1COOC°  K.  Yos  [1963] 
extended  the  computations  of  the  properties  of  hydx'ogen,  nitrogen,  oxygen 
and  air  into  the  ionization  regime.  Fay  [1962]  has  suggested  an  approxi¬ 
mate  mlxtui'e  rule  which  \m.s  used  by  Camac,  Fay,  Feinberg  erd  Kemp  [3  963] 
to  compute  the  properties  of  argon  in  the  ionized  region.  Other  results 
and  references  are  given  in  Pallone  and  van  Tassell  [1963]  and  Ahtye  [196^1]. 

A  defect  of  the  formulas  used  in  the  calculations  raer.tioned  e  jve  is 
t  lack  of  a  priori  agree.ment  of  the  thermal  conductivity  in  the  limit 
of  complete  ionization  wit'n  the  results  of  Spltzer  and  Harm  '1953]  oc 
landshoff  [1951]  for  the  fully  ionized  plasma.  This  defect  has  been 

These  two  books  will  hereafter  be  denoted  by  CC  and  MTGL 
respectively. 


realized  by  each  of  the  authors  and  they  have  modified  their  formulas  to 
obtain  pix>per  agreement  in  the  fully  ionized  limit.  Thus  Fay  [1962]  and 
Yos  [1963]  adjust  the  value  of  the  charged  particle  (Coulomb)  cross- 
section  in  their  approximate  fonaulas  to  obtain  agreement.  Eastlund 
[1963],  who  used  the  theoretical  formulas  of  MTGL,  inserts  some  constants 
to  insure  the  proper  I’miting  values.  Though  the  results  of  these  cal¬ 
culations  should  then  be  accurate  for  very  low  or  very  high  degrees  of 
ionization,  there  is  no  certainty  that  they  will  be  accurate  in  the 
intermediate  partially  ionized  "egion.  The  probable  accuracy  can  only 
be  determined  by  comparison  of  calculations  of  a  number  of  cafes  with 
resul  s  fr-om  rigorous  theoretical  formulas  which  "ive  the  correct  results 
at  low  and  high  ionization  without  any  ad  hoc  adjustment  of  constants. 

The  nature  of  this  difficulty  has  apparencly  also  been  realized  by 
Ahtye  [l96it-],  who  used  the  the  -etical  formulas  of  ITTGL  to  ccmput..  -he 
transport  coefficients  of  partially  ionized  argon.  We  will  see  later 
that  these  formulas,  which  are  really  Just  the  second  (and  first)  approxi¬ 
mation  to  the  exact  solutions,  need  to  be  extended  to  higher  approximations 
to  obtain  accii'’ate  values  of  the  thermal  conduct' vity  (and  viscosity)  for 
partially  ionized  gases. 


This  report  will  be  concerned  explicitly  '..’ith  the  calculation  of 
transrort  properties  of  partially  and  multiply  ionized  monatomic  gases, 
though  some  of  the  formulas  can  probably  be  used  with  -uccoss  for  poly¬ 
atomic  and  dissociit ; r.g  gases.  In  this  report  e>rrressjons  for  the  theimiai 
conduct i'/ity,  thermal  diffusion  and  the  multl-CGmpor.en*-  diff'usicn  coeffi¬ 
cients  of  monatomic  gas  mixtures  will  be  extenued  to  the  fourth  approxi¬ 
mation  in  the  Chapman- Ensh'^g  tncory.  The  visccsity  u^ll  be  extended  to 
-he  second  approximation  This  level  of  approximati''''n  has  been  considered 
previously  only  for  the  binary  gas  mixture  [Mason,  195Ta;  Saxena  and 
Joshi,  1963]  After  some  diocussion  of  the  charged  particle  cross- 
ssctions  these  exnrer^ions  will  be  applied  to  the  el ectron-singly  charged 
ion  plasma.  We  will  see  here  that  it  is  necessary'  to  use  at  least  the 
third  approximation  for  some  of  the  transport  properties  for  the  fully 
ionized  gas.  The  viscosity  also  is  In  error  by  at  least  15'^  if  only  the 
first  approximation  is  used.  At  this  point  we  can  examine  the  convergence 


of  the  approximations 
acMBe  simple  mixtures, 
to  the  calculation  of 
Ionized  argon  plasma. 


for  several  other  inverse-power  potentials  for 
Ir  1  later  p^per  these  expressions  will  be  applied 
the  transport  properties  of  singly  and  multiply 


2.  HIGHER  APPROXIMATIONS  IN  THE  CHAPMAN- ENSKOG  THEORY 


In  this  section  the  general  expressions  for  the  transport  properties 
for  aulti component  monatomic  gases  vill  be  extended  to  higher  approxima¬ 
tions.  Since  the  basic  theory'  is  already  well  covered  in  the  books  of 
Chapman  and  Cowling  [1952]  and  Hirschfelder,  Curtiss  and  Bird  [196U], 
only  a  cursory  review  will  be  given  in  this  section.  The  general  formu¬ 
lation  of  the  latter  reference  will  be  followed  in  this  report,  as  will 
the  majority  of  their  notation.  Th*’  numbers  in  brackets  will  refer  to 
r-juations  taken  directly  from  that  book. 

The  starting  point  for  the  Chapman- Enskog  theorj’'  is  the  multicompo¬ 
nent  Boltzmann  equation  for  the  distribution  function  of  the  ith  species 
[f^(r,v.,t)]  [7.1-25,10.1-5] 


df. 

1 

'ST 


j-1 


(2.1) 

The  right  side  of  this  equation  represents  the  rate  of  change  of  f^  due 
to  collisions  (or  encounters^  between  the  ith  species  and  itself  and  the 
other  v-1  species  of  the  gas.  It  has  been  written  here  in  two  equivalent 
forms.  Either  form  may  be  used  when  collisions  are  being  treated  classi¬ 
cally,  but  only  the  second  is  valid  when  quantum  mechanics  is  used.  The 
general  theory  follows  in  the  saime  way  for  either  case.  The  quantum 
mechanical  formulation  is  used  most  often  in  treating  electron-atom 
collisions  where  the  differencial  cross-section 
of  partial  wave  phase  shifts. 


A  list  of  symbols  is  given  on  page  vii. 
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is  obtained  in  terms 


Vs  note  that  this  equation  makes  no  provision  for  either  chemical 
reactions  or  excitation  of  internal  degrees  of  freedom  dxiring  a  collision. 
The  theory  is  then  strictly  applicable  to  non-reacting  monatomic  gases 
with  no  intemad  degrees  of  freedom.  The  rare  gases  at  temperatures 
below  about  5000*^K  and  above  those  temperatures  where  quantum  degeneracy 
is  important  come  closest  to  fulfilling  these  requirements.  When  ioniza¬ 
tion  takes  place  the  theory  is  no  longer  strictly  applicable  since  chemical 
reaction  and  excitation  of  electronic  states  are  important  processes  ta¬ 
king  place  during  collisions.  However,  only  a  very  small  percentage  of 
the  collisions  involve  excitation  or  chemical  reaction,  so  it  appears 
Justifiable  to  neglect  these  effects.  Except  in  the  case  of  the  thermal 
conductivity,  neglect  of  these  effects  has  been  successful  for  lower 
temperature  polyatomic  gases  where  rotational  and  vibrational  internal 
states  are  important.  Further  consideration  will  be  given  to  the  effect 
of  chemical  reaction  and  electronic  excitation  on  this  coefficient  in  a 
later  section. 


We  should  note  here  another  assumption  implicit  in  the  use  of  the 
Boltzmann  equation.  It  is  assumed,  namely',  that  it  is  possible  to  repre¬ 
sent  encounters  adequately  within  the  framework  of  this  equation.  It 
has  long  been  recognized  that  the  Boltzmann  collision  term  is  valid  for 
close  encounters  among  charged  particles,  i.e.  ft-r  those  for  which  the 
impact  parameter  b  is  of  the  order  of  the  mean  distance  of  closest 
approach 


2kT 


(2.2) 


here  written  for  the  encounter  of  two  particles  of  absolute  charge 

Z  e  and  Z  e  .  For  more  distant  encounters,  where  the  impact  parameter 
^  J 

is  of  the  order  of  or  greater  than  the  Debye  length 


the  Fc 

used  I 

derive 

Grad  [ 

collis 
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part  o 
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In  cas' 
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which  ( 
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with  c 
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This  definition  of  the  Debye  length  allows  for 
by  ions  as  well  as  by  electrons. 


plete  screening 
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the  Fokker-Planck  collision  operator  is  valid.  This  latter  form  has  beer, 
used  by  Cohen,  Spitzer  and  Routel^  [1950],  and  Spitzer  and  Hfirm  [19531  to 
derive  certain  of  the  transport  properties  of  fxilly- ionized  plasmas. 

Grad  [i960,  196la,  1961b]  ras  discussed  in  some  detail  the  problem  of  the 
collision  term  for  charged  particles  and  concluded  uhat,  at  least  to 
dominant  order,  i.e.  to  order  a  ,  where  A~  d/b^  ,  the  Fokker-Planck 
part  of  the  collision  terrr  may  be  replaced  by  the  Boltzmann  term  if  inte¬ 
grations  over  the  impact  parameter  b  are  cut  off  at  the  Debye  length. 

In  cases  where  ^  A  »  1  this  approach  is  satisfactory,  because  the 
neglected  terms  are  of  order  unity.  However,  a  large  number  of  plasmas 
which  are  presently  being  studied  do  not  satisfy  this  criterion.  For 
example,  cesium  vapor  at  1  mm  and  500o'^K  would  be  about  9C^^  ionized 

with  d~l.llxlO~^  cm  ,  n  ~  9.6xlO^Vcc  and  A~  5  ■  Argon  at  1  atm 

o  ^  —  6 

and  15000  K  has  about  oZ'*-  ionization  with  d~i.U2<10  cm  , 

n  9. 6x10 A~  36  •  Thus  we  see  that  teitas  of  order  unity 

should  really  be  included  in  considerations  of  the  charged  particle 

interactions  for  these  cases.  In  a  remarkable  paper  Kihara  and  Aono 

[1963]  have  accomplished  Just  this  step,  neglecting  only  terms  of  order 

A  or  smeJ-ler.  Their  method  appears  rather  difficult  to  apply  without 

going  into  considerable  denail,  but  we  can  use  an  approximate  treatment 

which,  for  the  cases  considered  by  Kihara,  Aono  and  Itikawa  [1963],  yivvs 

aLmost  the  same  results.  In  place  of  the  Coulomb  potential,  for  small 

angles  of  deflection  we  will  use  the  spherically  symmetric  shielded 

potential.  This  potential  has  been  considered  earlier  by  Liboff  [19591 

and  Ki.hara  [1959]-  Further  details  will  be  given  in  Section  3- 


2.1  Review  of  the  Char  man- Enskog  Method 

The  first  step  in  the  Chapman- Enskog  theory  is  to  assume  that  each 
distribution  function  can  he  written  in  the  form 


fi  =  fJ‘^^(l+Ti+  -) 


Jo] 


(2.I1) 


where  f^''"'  turns  out  to  he  the  Boltzmann  distribution  [?•  3-131 
/  m 


(2.5) 
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Here  and  are,  respectively,  the  number  density,  the  mass, 

and  the  velocity  of  the  ith  species,  is  the  mean  gas  velocity  and 

the  T  is  the  gas  temperature.  Upon  substitution  of  (2.5)  in  the  Boltz- 
mann  equation  (2.1),  it  is  found  that  the  perturbaoion  function  qp^  must 
be  of  the  form  (7.3-29,  35,36,3?] 


.  Aj(Wj)Wj-  ^  i  T)  :  ^  . 


V 

“I  (Wj) 


(2.6) 


where  d  will  be  defined  in  Section  2.2  . 

The  problem  of  determining  the  distribution  function,  and  hence  the 
transport  properties,  has  now  been  reduced  to  one  of  determining  the 
scalars  and  which  are  functions  of  the  magnitude  of  the 

reduced  velocity  ^i  ^i~'^i*^0  ‘ 

expand  each  of  these  scalar  functions  in  a  series  of  Sonine  (or  associ¬ 
ated  Laguerre)  polynomials  of  (7 •3-5'^].  The  and  functions 
are  expanded  in  terms  of  polynomials  of  order  3/2  and  the  B^  in 
polynomials  of  order  5/2  [7.3-6l], 


s-j. 

w  -T  V  ^i/2  '"f) 


5-J. 

®i^^’^i^  ^Im  ^/2 


(2.7) 


^  ^/2  (^i^ 


The  number  of  terms  retained  in  these  series  (i.e.  5-1)  will  be  called 

*) 

the  level  of  approximation.  Usually  one  or  two  terms  are  sufficient 


'This  convention  agrees  with  Chapter  ?  of  MTGL  but  disagrees  with 
Chapter  8  of  that  bo  k. 
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to  give  accurate  results  for  the  transport  coefficients.  This  is 
fortunate  because  the  formulas  beccrne  quite  cc.Tiplicated  at  only  the 
second  approxiniation.  We  vill  see  later  that  it  is  necessary  to  use  at 
least  the  third  apprcxirriaticn  when  computing  the  thermal  conductivity  and 
thermal  diffusion  coefficient  of  partially  ionized  gases,  though  the 
second  approximation  appears  to  be  quite  adequate  for  the  other  properties 
We  vill  now  proceed  to  develop  the  formula.*;  for  the  transport  coefficients 
Again  free  use  will  be  made  of  equations  already  given  in  MTGL. 


2.2  Diffusion  Coefficients 

The  general  expression  for  the  diffusion  velocity  of  the  ith 
species  relative  to  the  mass  average  velocity  is  given  by  [7-^-3] 


2  r- 

n  \ 


<  V  >  =  —  )  m  B.  .d  -  dJ 
1  n^P  ^  J  J  i 


(2.8) 


and  are  the  multicomponent  diffusion  and  the  thennal 

diffusion  coefficients,  respectively,  d  is  the  sum  of  the  concentra- 

J 

ticn,  pressure  and  external  driving  forces  for  diffusion  [7*3-27], 


J  d?\"  /  I"  =  /  a?  \  Pf-  /  "j  J  .4  ‘  ■ 


(2.9) 


where 


is  the  external  force  acting  on  the  £th  species  of  the  gas. 


In  terras  of  the  coefficients  of  expansion  of  the  distribution  function 
(2.7),  the  diffusion  coefficients  are  given  by  [7 *^-8,9] 


^  ,'2M  Ji 


ij  2nmj\^  i. 


(2.10) 


_T  ^i°^l  /2kf" 
i  ~  2  ^iC 


(2.11) 


We  can  now  examine  in  some  detail  the  solution  of  the  equations  for 
tlie  higher  approximations  to  the  multicomponent  diffusion  coefficient. 


The  fbnoula  for  this  coefficient  to  the  first  approxiniation  has  already 
been  given  in  MTGL,  and  the  same  general  procedxire  will  be  followed  as 

was  \xsed  there. 

hk 

The  coefficients  Cj^^  are  determined  from  the  set  of  equations 

[7.3-75] 


V 

IE 

J-1 


«nC  hk  ^  „ml  ^hk  ^  ^m2  hk  ^  „m3  ^hk 
"jO  "^ij  ‘^Jl  ^  "J2  "  ^ij  =J3 


1  „hk  _  3  IM 

J  ■  im  “  2\/ 


(S,,  -5..  )5  „ 
'  ik  ih  mC 


(2.1?) 

where  m=sO, 1,2,3  (for  the  fourth  approximation)  and  i=l,2,3‘“V  .  The 

solution  of  these  equations  may  be  obtained  with  the  aid  of  Cramer' s 

rule  (Hildebrand,  1952]  as  the  ratio  of  two  2v  x  2v  determinants. 

The  denominator  is  the  determinant  of  the  coefficients  Q™  ,  and  the 

J 

numerator  involves  the  coefficients  with  one  coluian  replaced  by  the 
right-hand  side  of  Eq.  (2.12).  A  more  convenient  form,  equivalent  to 
that  obtained  with  Cramer's  rule,  has  been  introduced  in  MTGL.  Before 
proceeding  to  the  explicit  expression,  it  is  convenient  to  define  new 
elements  different  from  those  used  in  MTGL.  In  this  way  the  number  of 
operations  performed  in  the  calculation  of  the  various  coefficients  will 
be  materially  reduced.  We  let 


..  FT 


Siim-  Kir 

“wT  ^iO  ’^^^ih"^lk^ 


(2.13) 


and  the  solution  for  the  fourth  approximation  to  the  multicomponent 
diffusion  coefficients  is  then 
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where  |ql  is  the  deteminant  forced  from  the  ri’umerstor  by  deleting 

the  last  row  and  last  column.  (For  convenience  of  notation  the  s^-mbol 
mm' 

q^^j  represents  a  block  of  elements  v;ith  both  i  and  j  ranging  from 
1  to  V  .)  To  obtain  the  first,  second  and  third  approximations  to 
we  simply  delete  all  q  "blocks"  except  those  with  m=n'=05m, 
m'^  ljni,m'5  2  respectively,  from  both  the  numerator  and  denominator. 

mm' 

iiefore  proceeding  to  explicit  formulas  for  the  elements,  we 

can  easily  derive  the  form  of  the  thermal  diffusion  coefficient.  For  the 
fourth  approximation  the  a  are  solutions  of  the  equations  [7.3-75] 

1(J 


it 


^  ^ml  ^  _m2  ^  .^m3  □ 

+  ft-  ^^.-1  ft-  ft.  .^1-.  =  -  Rj 

^ij  jO  jl  J2  J3j  i 


^5°i  /2kT 


where  i-1,2, -v  ,  and  m=0,l,2,3  •  Here  R.  =  —i —  / -  6.  .  We 

1.II1  ^  ^  Xm 


can  solve  for  a^^  as  in  the  previous  case  and  write  the  thermal  diffu¬ 
sion  coefficient  as 
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,-T,  15 

■  TiT  — fqj — 


(2.15) 


To  obtain  the  third  and  second  approximations  to  the  thermal  diffusion 

rnni* 

coefficient  we  simply  delete  all  blocks  containing  with,  respec¬ 

tively,  m  or  m'=3  and  m  or  m' i  2  .  In  the  first  approximation 
T 

Is  Identically  zero. 

The  general  expression  for  the  determinant  elements  occurring  in 
the  above  relations  is  f7*3-7l] 


r  “il 


Here  the  are  the  Sonine  polynomials  of  mentioned  earlier 

and  defined  in  [7.3-57]*  The  expressions  in  brackets  are  defined  in 
[7*3-^3]*  These  "bracket  integrals"  can  be  reduced  to  a  simpler  form  in 
tenns  of  the  Integrals  [8.2-8] 


^  2  (i>s)*^  - kjj^ 
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(3+l)l[2i+l-(-l)^] 
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qH  (g)  =  2ir  \  o,  ,(x)(l-  cos'®  x)  sin  X  <3x  =  27r  \  (l-  cos'®  x)'^^ 


(2.18) 
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Here  a.  .(x)  is  the  differential  cross-section,  x  is  the  angle  of 

2  2 

deflection  during  a  collision,  and  y  =l/2  ((^ .  .g  /kT)  ,  y  being  the 

J  ^  J 

reduced  mass  of  the  two  species  i  and  j  . 

General  formulas  for  this  reduction  have  been  given  in  CC.  Some  of 
the  explicit  formulas  are  given  in  CC  and  NfTGL.  Other^^  have  Cden  derived 
from  the  work,  of  M^ason  [1957a.].  Relations  for  obtaining  the  bracket 
expressions  from  his  work  are  collected  in  the  appendix.  With  the  aid 
of  these  formulas  and  af^^r  some  algebraic  manipulation,  the  following 
expressions  for  the  q™  elements  may  be  derived; 

V  l/  2  r*  ■ 

<>00  00  ^00  oV  '^^“i  _2  .(I,!)""  ri  e  \ 
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(2.19c) 
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2.3  Energy  Flow  and  Thermal  Conductivity 


The  genei*al  expression  for  the  energy  flux  in  a  monatomic,  non- 
reacting  gas  with  rc  internal  degrees  of  freedom  is  given  by  [7.^-3C] 


V  V 

q  =  |kTy  n  <  V  >  -  X'  —  -  nkT  V  d 

2  Z.  J  J  /_  n  m  j  J 

J=1  J=1  ^  ^ 


(2.20) 


The  first  term  in  this  expression  represents  the  translational  enthalpy 

carried  by  each  species  as  it  diffuses  relative  to  the  mass-average 

velocity  of  the  gas  mixture.  The  last  term  represents  a  coupling  between 

diffusion  forces  and  heat  flow.  It  is  conventional  to  eliminate  the 

explicit  occurrence  of  d  in  this  expression  by  using  the  formula  for 

J 

the  diffusion  velocity  (2.3)  which  can  be  rewritten  in  the  following  form 


i  W:  “  - 


°lC  <  '^1  ^  » 

2  2  m 

n  n  m^T  or 


(Dj^=0)  (2.21) 


In  order  to  eliminate  the  d  from  (2.20)  we  can  consider  (2.21)  as  a 

set  of  linear  equations  to  be  solved  for  the  uPikncwns  d^  .  Defining 

E  as  an  element  of  the  inverse  of  the  matrix  whose  general  element  is 
J  ^ 

m.D.  ve  can  wr^te  a  formal  solution  for  d  as 

J  i  J  J 


(2.22) 


Then  we  have 
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The  energy  flux  then  becomes 


V  T 
E. .D: 


q  =  |kTy  l-IS-V  JLli-  n  <  V  >  -  X  — 

•*8  L  5  n  J  J 


where  the  coefficient  cf  thermal  conductivity  is  given  by 


(2.23) 


X  *  X’ 


Vi 


E  0*^0*^ 
i  J 


(2.24) 


Since  no  chemical  reactions  or  internal  transport  of  energy  are  considered 
’in  this  expression,  this  X  is  often  called  the  frozen  thermal  conduc¬ 
tivity  and  is  sometimes  denoted  by  X^  . 

The  solution  above  for  d  has  been  called  a  formal  solution  since 

V 

no  explicit  formula  has  been  given  for  finding  the  inverse  matrix  £  in 

inni' 

terms  of  the  elements  q.  .  In  actual  calculations  the  diffusion  coef- 

Ij 

flcients  D.  .  will  be  calculated  and  then  £  =  (Dm)  [written  here 
in  matrix  form]  may  be  found  from  a  standard  algorittmi.  If  the  calcu¬ 
lations  are  dene  on  a  computer  such  as  the  IBM  7090 >  which  is  used  for 
such  calculations  at  Stanford,  then  there  is  generally  a  library  proce¬ 
dure  available  which  can  be  used  to  accomplish  the  inversion. 

T 

Expressions  for  D  and  D  were  given  in  the  last  section.  X' 

*•  1 J 

is  given  in  terms  of  the  expansion  coefficients  of  the  scalar  function 
A^(W^)  by  [7.i^-35] 


X'  =  - 
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J  \  *“4 

.i»l 


(2.25) 


where  a^^  is  obtained  f 'om  the  equations  [7*3-33]  already  given  in 
connection  with  the  thermal  diffusion  coefficients.  We  can  solve  for 
the  ajj_  Just  as  in  that  section,  X'  then  being  given  by  a  sum  of 
such  expressions.  The  final  result  is 
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30  1 

1 

1 

1 

rf33  1 

I 

0 

0  -1 

\ 

0  1 

1 

0  1 

0 

i 

1 

This  is  the  f :rEula  for  X'  in  the  fourth  approximation,  anJ  to 

obtain  the  third  :r  second  approximation  we  simply  delete  the  blocks 
nun* 

involving  q.  .  vii.h,  respectively,  m  or  m’=3  and  m  or  m'  >  2  . 

J 

This  is  the  same  7rescrirwion  as  given  in  the  section  on  the  diffusion 

coefficients.  Trouble  arises  when  we  attempt  to  obtain  a  first  approxi- 

nun* 

mation  by  deleting  also  the  q. .  's  Involving  m  or  m' =1  ,  for  then 

^  J 

the  determinant  in  the  r/umerator  vanishes.  This  is  equivalent  to 
retaining  only  the  first  term  in  the  expansion  of  A^(W^)  .  We  can 
define  a  first  approximation  to  the  thermal  conductivity  as  that  obtained 
when  we  keep  only  the  second  term  in  the  expansion  of  A^(W^).  [See 
Eq.  (2.7).]  The  formula  for  the  first  approximation  can  then  be  obtained 
from  (2.26)  by  deleting  all  q  blocks  except  those  with  m=m'=l  . 

Muckenfuss  and  Curtiss  [1958]  have  shewn  that  the  first  approxi¬ 
mation  obtained  in  this  way  is  identical  with  the  total  thermal  conduc¬ 
tivity  to  the  second  apprcxinaticn.  This  arises  in  part  because  in  the 

T 

first  approximation  [X]^=[X’]^  •  From  irreversible  thermo¬ 

dynamics  it  is  always  true  that  X  <  X'  and  they  show  that  the  difference 
is  sufficient  to  make  =  [X^^-  The  complete  formula  to  the  second 

approximation  is  then 
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(2.27) 


ny/m^  0 


ThlB  form  is  considerably  simpler  than  the  one  given  earlier  for  the 
second  approx irsat ion  and  is  usually  used  in  "rigorous"  calculations  of 
the  thermal  conductivity.  We  will  see  later  that  this  approximation  gives 
answrs  which  are  much  too  low  for  the  fully  ionized  gas,  and  presumably 
also  for  the  partially  ionized  gas.  It  is  necessary  then  to  go  to  the 
third  approximation  for  accurate  values  of  this  coefficient.  Unfortu¬ 
nately  no  reduction  similar  to  the  above  exists  at  present  for  the  third 
or  higher  approximations. 

The  fora  (2.23)  of  the  heat  flux  is  most  convenient  when  the  thermal 
conductivity  is  defined  as  the  coefficient  of  the  temperature  gradient  in 
the  expression  for  the  flux  when  no  diffusion  is  taking  place.  In  (2.23) 
we  recognize  X  as  the  true  thermal  conductivity.  However,  (2.23)  is 
resLlly  a  mixed  type  of  expression  in  that  the  heat  flux  is  given  in  terms 
of  a  force,  the  temperature  gradient,  and  other  fluxes.  If  we  wish  to 
compute  the  total  heat  flux  at  a  point  in  the  gas  and  are  not  concerned 
with  the  "true"  thermal  conductivity,  then  this  expression  is  often 
not-  the  most  convenient  fora.  It  is  often  more  practical  to  write  che 
expression  for  the  heat  flux  from  the  outset  solely  in  terms  of  the 

)  -N— ,  i 

forces  d.  and  dT/or  .  This  is  the  convention  followed  in  irre- 
J 

versible  thermodynamics. 

To  write  the  heat  flux  in  terms  of  forces  we  use  Eq.  (2.8)  to 
eliminate  the  number  fluxes  from  (2.20).  After  some  manipulation  we 


arrive  at 
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(2.28) 


actually  has  the  dimensions  of  L  but  it  is  convenient  to 


refer  to  it  as  a  force. 
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is  the  total  enthalpy  of  the  ith  species  (see  below).  The 


where  h^ 

disadvantage  of  this  form  nay  be  seen  if  we  try  to  find  the  usual  coef¬ 
ficient  of  thermal  conductivity.  To  do  so  we  must  set  the  diffusion 
fluxes  equal  to  zero  in  (2.8),  solve  for  the  diffusion  forces  d.  ,  and 

V 

eliminate  the  d.  in  (2.28)  with  the  resultant  expressions.  Clearly 
(2.23)  and  (2.24)  are  more  convenient  forms  for  this  purpose,  even 
though  they  require  finding  the  inverse  of  Dm. 


The  formulas  above  for  the  energy  or  heat  flux  vector  apply  still 
only  to  a  non-reacting,  monatomic  gas  with  no  internal  degrees  of 
freedom.  Neglecting  the  effect  of  the  chemical  reactions  on  the  distri¬ 
bution  functions  f  ,  their  only  effect  is  to  add  an  energy  of  formation 

1^.  i 

h^  to  the  translational  enthalpy  carried  by  each  species  as  it  diffuses 
relative  to  the  mean  gas  motion.  Thus  we  replace  5/2  kT  in  (2.23)  by 
[  (5/2)  (kT/m  )+h*^Jm  and  bring  it  inside  the  summation  over  J  , 

J  J  J 


(2.29) 


In  this  expression  the  second  term  multiplying  the  diffusion  velocities 
in  (2.23)  has  been  dropped.  Calculations  show  this  term  to  be  negligible 
compared  to  the  first  term. 

Butler  and  Brokav  [1957;  also,  Brokaw,  196c]  have  carried  out  a 
reduction  of  (2.29)  in  the  case  of  chemical  equilibrium  to  give  l.  much 
simpler  form 


q  = 


react 


(2.30) 


To  derive  this  relation  they  neglect  thermal  diffusion  and  use  only  the 
first  approximation  to  the  diffusion  coefficients.  We  will  see  later 
that  the  latter  can  be  in  error  up  to  aliaost  a  factor  of  2  in  the  ionized 
gas.  Also  the  thermal  diffusion  does  play  an  important  role  when  free 
electrons  are  present  and  -^rrbably  cannot  be  neglected  in  this  case.  It 
thus  appears  that,  even  if  chemical  eq-’lllbriun  did  exist,  the  formula 
which  they  derive  would  not  be  applicable  to  the  partially  ionized  gas. 
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It  seems  best,  then,  to  leave  the  expression  for  the  heat  flow  in  the  form 
of  (2.23)>  which  is  also  more  general  than  the  expression  they  derive 
since  it  applies  eq.ually  to  situations  with  chemical  non-equilibrium. 

The  monatomic  gases  are  generally  considered  to  be  without  internal 
degrees  of  freedom  for  purposes  of  calculations  of  thermal  conductivity. 
Fbr  argon,  for  example,  at  temperatures  below  about  15C00°K  this  turns 
out  to  be  a  very  good  approximation.  Results  of  Drellishak,  Knopp, 
and  Cambel  [1963I  show  that  at  about  this  temperature  the  population  of 
the  excited  states  is  sxifficient  to  cause  the  atom  partition  function  to 
begin  to  deviate  from  the  unexcited  value  of  1  .  Similar  deviations 
would  be  expected  for  alkali  metal  vapors  at  a  much  lower  temperature. 
When  electronic  excitation  becomes  important,  the  possibility  of  heat 
transport  by  excited  states  must  be  considered.  A  similar  problem  arises 
in  polyatomic  gases  at  lower  temperatures  where  energy  transport  by  ro¬ 
tational  and  vibrational  states  must  be  considered.  This  latter  case 
has  been  handled  with  moderate  success  by  the  Eucken  factor,  which  is 
an  approximate  correction  to  account  for  this  form  of  energy  transport. 

The  first  step  in  allowing  for  excitation  of  electronic  states  is 
to  r')write  (2.29)  as 


■I" 

J=1 


V.  ^  ,  3T 

.m.h,  <  V.>  -  X  — 

J  J  J  )  g; 


(2.29^ 


where  the  total  enthalpy  of  each  species  h  now  includes  the  transla 


tional  enthalpy  {5/2)(k/ra  )T  ,  the  energy  of  formation  h^  ,  and  the 
energy  stored  in  excited  electronic  states.  In  practice,  the  energy  of 
formation  (reactive  energy)  will  probably  swamp  the  other  two  forms. 


Now  because  excited  species  can  diffuse  relative  to  the  unexcited 
species,  thus  transporting  energy  of  ■'  rcitation,  a  correction  must 
generally  be  applied  to  the  coefficient  of  thermal  conductivity.  For 
polyatomic  gases  this  correction  takos  the  form  of  an  Eucken- type 
correction  [Mason  and  tonchick,  19^2;  Hirschfelder,  1957a,b].  How¬ 
ever,  Hirschfelder  has  pointed  out  that  the  diffusion  coefficient  for 
an  electronically  excited  atom  through  atoms  in  the  ground  state  is 
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probably  very  small.  He  arrives  at  this  conclusion  from  estimates  of 
the  collision  diameter  for  excited  species.  For  example,  argon  in  the 
first  excited  state  would  have  a  collision  diameter  about  3  times  that 
of  the  ground  state.  Mason,  Vanderslice  and  Yos  [19591  have  pointed  out 
that  the  diffusion  cross-section  would  be  further  reduced  because  of 
resonance  exchange  collisions  which  can  occur  with  atoms  in  the  ground 
state.  It  seems  likely,  then,  that  we  can  neglect  altogether  the  effect 
of  excited  electronic  states  on  the  coefficient  of  thermal  conductivity 
of  monatomic  gases. 


2.1-  Viscosity 

In  terms  of  the  coefficients  of  expansion  of  E^(W^)  [see  Eq.  (1.7) ]> 
the  viscosity  Is  given  by  [7-^-2^- 1 
V 

J=1 

Since  the  convergence  of  the  series  for  the  viscosity  is  typically  quite 
rapid,  we  will  consider  only  the  second  apprcxi.mation.  The  equations  to 
solve  for  b,,^  are  [7 •5-76] 


V 


ij  JO 


''ll 
i  J  Jl 


where  i=l,2,-v  . 


(2.32) 


Tne  solution  to  these  equations  can  be  obtained  just  as  it  was  for 
X'  occurring  in  the  theicnal  conductivity.  As  before  we  define  new 
determinant  elements  to  simplify  calculations 


"'mm'  1^1 

iij  's'-w  '*ij 


(2.33) 
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aud  arrive  at  the  formula 


hi 


5  Vamd* 
2  “  “  2  “RT 


'XtO  I  "01  I  r 


"10  I  "11  I  ^ 


I  0  I  0 


(2.34) 


To  obtain  the  first  approximation  from  this  expression  we  delete  the 
blocks  with  m  cr  m'=l  .  The  general  expression  for  the  determinant 
elements  is  (7.3-71] 


“iKjLV5/2  ^l®5/2 

1-1  I 


"  »J/  Pl^/2  >  ^5/2 


u 


J 


with 


(2.35) 


Wi  =  ■  3  ^  • 


(2.36) 


As  for  the  determinant  elements  occurring  in  the  expressions  for  the 
other  transport  coefficients,  these  may  be  written  entirely  in  terms  of 
the  integrals.  Making  use  of  the  formulas  given  in  NfTGL  and 

in  Saxena  and  Joshi  [1963]  ve  find 


2^ 


1/2 

I 


/^\y 


3/2  ^  "ii  -  “1/ 


ij  jr  j 


^  "u 


(2.37a) 
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v2  V  3/2  r 

«1J  -  l^jZ  (  )5/2  ™u 

X*“-i  1  Xf 


;  B  «■“•  -  ";:■”•)•  ■,  »,rt. >('«■'’•■  “S’”')} 

(2.37b) 


"10  r^Voi 

"ij  =  [iK 


(2.37c) 


/  V  2  V  1/2  r 

"11  q  \\  V  ^i”.£  _  2  L  s 

•  j  "j  ('■'"“j  *  2‘‘5=>!)  “u  * 

*  J  ^(15^  +  WmJ)  n'2'®>*  -  56m3  n<2’3)*  , 

+  k)in^  I 

i  ii 

-j 


(2.37d) 


2.5  Electrical  Properties 

The  equations  and  properties  ^iven  in  the  previous  sections  are 

sufficient  to  completely  specify  the  mass  (or  number j,  mor'entum,  and 

energy  fluxes  in  a  mixture  of  monatomic  gases.  T'c  is  worthwhile  to 

examine  the  special  form  which  these  equations  take  when  electrons  and 

*) 

ions  are  pre-ent  and  when  there  is  an  external  electric  field.  In 

this  way  we  will  derive  an  explicit  expression  for  the  electrical  con¬ 
ductivity.  Without  loss  of  generality  we  can  consider  for  the  moment  a 
mixture  without  concentration  or  pressure  gradients.  Equation  (2.9)  for 


No  magnetic  fields  are  considered  in  the  present  treatment. 
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the  dlffusiori  forces  becomes 


a.m, 

J  J  P 


d  E.J  -V  .x 

j  nkTp  m 

^  k=l 


(2.38) 


Let  species  1  be  the  electron,  species  2  through  ^  be  ions,  each 
vith  charge  eZ^  ,  and  the  remaining  v-^  particles  be  neutral.  Then 
the  external  forces  are 


=-eE 


=  eZ^E 


2  <  1  <  ? 


(2.39) 


='l'° 


1  >  C 


A  \\  *  •  =  0 


since,  for  a  neutral  plasma,  n., 


Z  .  The  diff’jsion  forces  are 


-•  n-j^eE 
*^1  ”  nkT 


n,Z  eE 

X  1 


i  >  ^ 


2  <  i  <  C 


(2.40) 


Using  these  expressions  in  Eq.  (2.8)  the  electro.,  and  ion  contributions 
to  the  current  are 


d  is  really  not  a  force.  See  footnote  on  page  19 . 
J 


ne^ 


1 


en,  <  V.>  Z. 
i  i  1 


"A'li 


n.m.zp.  , 
J  J  J  ij 


i2.kl) 
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The  electrical  conductivity  is  defined  as  \'ie  coefficient  of  E 
in  the  expression  for  the  total  current, 


®  *  in?r 


n.n.Z.D^^  4  riTiiip.D^  "  n.m.Z.D,. 

JJJIJ  IIJJI  iiiJi 

i=2 


(2.42) 


Nov  m.  «  m  so  we  can  examine  this  expression  for  terms  of  various 

orders  of  m.  t  .  The  diffusion  coefficients  P  are  of  order 
J  ^  J 

(u,. so  (l/m^^)  ,  (l/m^'^^)  .  Thus  the 

1/2 

first  term  in  rhis  expression  is  of  order  >  the  second  of 


1/2  1/2 

^rder  and  the  third  of  order  2p 

1  U  J 


J'"l  '  . 

Except  in  cases  of  large 


Z,  ,  only  the  first  term  is  important  and  we  make  the  approximation 


2 

e  n 
^  ~  pkT 


(2.42a) 


The  coefficient  of  the  temperature  gradient  in  (2.41)  is  also  of 


interest.  The  complete  expression  is 


T  m,  ^ 

J=2 


and  to  highest  order  in  the  masses, 
T 

“  Tm^ 


(2.43) 


(2.43a) 
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Again  we  examine  the  order  of  the  various  terms  in  these  expressions. 

In  (2.45)  ve  would  be  tempted  to  say  that  and  h^(l/m^)  so 

that  all  of  these  terns  except  for  the  one  corresponding  to  the  electrons 
can  be  neglected.  However,  if  chemical  reactions  are  possible  -  and 
ionization  or  recombination  are  always  possible  when  electrons  are 
present  -  then  h^  might  also  include  a  very  large  ionization  energy. 
Only  In  the  case  of  a  flow  situation  so  rapid  that  the  chemical  compo¬ 
sition  is  essentially  frozen  would  this  type  of  reduction  by  order  in 
mass  be  permissible  without  more  careful  examination  of  each  h^  . 
this  reservation  in  mind  we  go  ed'ead  and  make  the  reduction  to  get 


With 


The  reductiori  of  (2.46)  follows  in  the  same  aianaer  with  the  above  noted 

,  the  second  of  order 


1/2 


reservation.  The  first  term  is  of  order 

1'2  1/2  ^ 

^7  >  and  the  fourth  of  order  a/"  .  In  the  terms  involving  the  ther- 

-i- 

aal  diffusion  coefficient,  the  electron  term  is  still  dominant.  The 
simplified  form  of  (2.U6)  is 


0  = 


3  en  ^ 

2  P  _ 

J=2 


ec: 


(2.46a) 


The  expressions  for  the  heat  and  current  flow  in  the  presence  of 
zhe  electric  field  are  new 
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hr 
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K-  ^ 
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(2.47) 


To  get  the  true  thermal  conductivity  K 
tc  get 


K  =  €K' 


€  =  1 


.  ^ 

CK' 


we  set  J=0  and  solve  for  E 


(2.48) 


The  modification  of  these  expressions  to  include  concentration  and 
pressure  diffusion  forces  follows  in  a  straightforward  manner  from 
(2.9)*  The  complete  formulas  will  not  be  given  here. 
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3.  APPLICATIONS 


We  are  now  in  a  position  to  compute  the  transport  properties  to  the 
fourth  approximation  for  all  except  the  viscosity,  and  this  to  ohe  seccni 
approximation.  We  note  that  the  complexity  of  the  formulas  is  increased 
manifold  with  each  increase  in  level  of  anvroximaticn.  Further,  each 

(i ) 

added  level  requires  the  computation  of  a  kinetic  cross-section  Q 

2  K  Ji  S  )  * 

for  a  new  £  ,  and  two  new  integrals,  v  a  Cl  ,  for  each  lower 

value  of  £  .  Before  attempting  to  calculate  the  properties  of  a  real 
mixture  then,  it  is  important  to  have  an  idea  of  the  level  of  approxi¬ 
mation  necessary  to  obtain  reasonably  accurate  results.  For  the  lavs 
of  force  effective  in  low  temperature  gas  mixtures,  the  second  approxi¬ 
mation  for  the  thermal  conductivity  is  quite  adequate  and  the  thermal 
diffusion  is  usually  ignored.  For  the  ordinary'  diffusion  coefficient 
and  the  viscosity,  the  first  approximation  is  accurate  enough.  These  are 
Just  the  levels  of  approximation  for  which  expressions  are  given  In  IfTlL. 

In  the  case  of  a  partially  ionized  plasma  there  are  several  forces 
operative  between  the  various  pairs  of  particles,  and  it  is  difficult 
to  examine  the  convergence  fcr  this  mixture  without  some  particular  gas 
in  mind.  We  can,  however,  examine  the  convergence  of  the  limiting  case 
of  a  fully-ionized  binary  gas,  since  only  one  force  law  is  needed  here. 
This  special  case  has  already  been  treated  with  various  simplifying 
assumptions  by  several  other  workers,  and  we  can  compare  their  results 
with  those  obtained  with  the  expressions  of  this  report .  Before  pro¬ 
ceeding  to  this  example,  we  should  examine  the  charged  particle  inter¬ 
action  with  the  idea  of  re'aining  higher  order  terms  in  the  cross- 
section.  In  the  first  section, expressions  for  the  cross-sections 
necessary  for  up  to  the  fo’urth  approximation  will  be  developed  following 
the  method  used  by  Liboff  [1/59]-  In  the  pi*esent  report  only  the  domi¬ 
nant  log  term  in  these  cross-sections  will  be  used,  since  the  effect  of 
the  higher  order  term  of  order  unity  depends  on  the  temperature  and 
charged  particle  density.  Hence  its  effect  on  the  transport  properties 
must  be  Judged  in  particular  cases. 
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Two  methods  for  calculating  the  properties  of  the  binary  plasma 
suggest  themselves.  One  would  be  to  consider  the  special  form  that  the 
general  expressions  take  in  this  case,  and  to  neglect  terms  of  low  order 
in  the  mass  ratio  to  simplify  the  calculations.  Landshoff  [19^9>  1951] 
has  essentially  used  this  method,  starting  from  the  somewhat  different 
but  equivalent  formulation  of  CC  for  the  binary  gas.  In  the  present 
case  it  was  decided  not  to  consider  the  reductions  possible  in  special 
cases.  Rather,  the  formulas  were  programmed  for  a  computer  in  their  most 
general  form  and  special  cases  were  computed  with  this  general  program. 
The  program  has  been  written  in  Subalgol,  a  version  of  Balgol  developed 
at  Stanford  for  the  IKM  7090-  As  input  to  the  program  it  is  necessary 
to  specify  the  masses  and  number  densities  of  each  species,  and  pro¬ 
cedures  for  computing  the  various  cross-sections. 

After  the  results  for  the  binary  plasma  have  been  compared  with 
those  from  other  sources,  convergence  will  be  examined  for  three  ether 

special  mixtures.  Consideration  of  these  cases  is  fruitful  in  two  ways, 
‘•lest  important,  it  demonstrates  the  range  of  rates  of  convergence  which 
can  be  expected  for  real  gas  mixtures.  Also,  it  c ffers  a  check  on  the 
expressions  of  Section  1  and  their  translation  to  the  computer  program. 

In  some  of  the  cases  to  be  considered  the  results  may  be  compared  with 
the  exact  answer;  in  other  cases  the  rate  of  convergence  has  been 
derived  by  different  methods. 


3. 1  Charged  Particle  Cross-Sections 


The  usual  procedure  for  the  charged  particle  cross-section  is  to 
take  the  Coulomb  law  of  interact* on 


(3.1) 


as  valid  up  to  a  finite  distance,  either  the  Debye  length  d  [Grad, 
1961b]  or  the  laterelectron  distance  h^n^^^^  [CC,  p.  l'^9]*  and  to 
neglect  all  interactions  beyond  this  distance.  From  the  discussion  of 
Section  II,  It  is  '■oparent  that  for  a  large  variety  of  plasmas  the  next 
higher  order  terra  beyond  the  dominant  x/b^  where  x:=d  or  h  , 
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should  be  considered.  This  could  be  accomplishod  v^th  the  usual  cut-off 
procedure,  but  a  more  realistic  method  wo\ild  be  to  use  the  Coulomb  poten¬ 
tial  for  close  encounters  and  a  shielded  Coulomb  potential 


Z  Z  e 

i  J 


(3.2) 


for  the  more  distant  encounters.  Liboff  [1959)  has  evaluated  some  of  the 
necessary  integrals  with  this  method.  We  can  follow  his  method  to  cc^m- 
pute  the  additional  integrals  necessary  for  the  level  of  approximation 
considered  here.  Later  some  of  the  results  ’.d.il  oe  compared  with  the 
more  exact  results  of  Kihara,  Aono  and  Itikava  [1963]. 

The  cross-section  integral  of  (2.18)  is  separated  into  two  parts 

Q^j^g)  =  +  27r^(i-  cos^  x)bdb  , 

the  first  to  be  computed  using  the  Coulomb  potential  (3*1)  and  the 
second  with  the  shielded  potential  (3.2).  L  is  some  distance  satisfying 
b^  «  L  «  d  where  b^  and  d  are  defined  in  (2.2)  and  (2.3)*  For 
the  Coulomb  potential  the  angle  of  deflection  is  related  to  the  impact 
parameter  b  by 


x^-1 

cos  X  =  -p— 
X  +1 


(3.3) 


where 


X  =  2b 


2 

Z  Z 

if 


Substituting  this  form  into  the  cross-section  integral  for  the  cases 
i«3,4  we  get  ^ 

..(3)  37r  ,_.x2  i. ,,,  ..  2x  2 


q(3)  .  2E  ^^,  2^  ^  2 

^  ^  °  ^  (x^l)^ 
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vith 


3 

0  2kT 


the  croi?s- sections  are  (including  those  obtained  by  Llboff) 


=  Utt  j  — j  en  j  -  1  ^  2 

/  Vq  2 


=  Stt  (^  )  &n  j  -  1  -  >  + 


=  127r[~J  ^(?^)-r->+^2 

/  voj  ° 
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(3.9) 


2  (  £  s )  * 

We  now  need  the  average  cross-sections  t  a  H'  ^  given  by  (2.17), 


which  can  be  rewritten  in  the  form 


00 

^  e'^  y®‘"^Q^^^y)dy  (3.10) 


Inserting  the  expressions  for  Q'  '(y)  we  arrive  at  integrals  of  the 
form  (csrconstant) 


r  ^  ^b^j  ■  =  r'(s)jt(s)  -  c  (3.11) 

s-1 

where  for  s  integer  the  Psi-function  ^(5)=  -7+  y  -  where  >=0.5772- 

Z_  n 

n=l 

is  Euler's  constant.  This  first  integral  may  be  evaluated  by  integration 
by  parts  or  located  in  Bierens  de  Haan  [1957].  The  general  expressions 


for  the  average  cross-sections  are 

r 


l,2U  A  -  i  -  2>  .  y  (i 


(3.12) 
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We  now  have  the  cross-sections  necessary  for  the  calculation  of  the 

transport  properties  of  mixtures  of  charged  particles  to  the  levels  of 

approximation  considered  in  Section  2.  To  dominant  order,  i.e.  when 

in  A  »  1  ,  these  results,  using  a  screened  Coulomb  potential,  agree  with 

those  obtained  by  using  the  unscreened  potential  together  with  a  cut-off 

of  the  integrals  at  the  Debye  length.  We  can  now  compare  some  of  these 

expressions  with  those  obtained  by  the  more  exact  method  of  Kihara, 

Aono  and  Itltawa  [1963].  They  examine  in  detail  several  phenomena  in  a 

binary  plasma  and  obtain  expressions  equivalent  to  the  average  cross- 
2  { £  s  ^  ^ 

sections  to  with  i=s=l  or  2  .  For  the  case  cf  relaxation 

between  electron  and  ion  temperatures  (£=s=l)  they  obtain  the  same 
expression  as  in  (3*12)  with  a  different  screening  constant 


j2  kT 
®  47m 

e 


(3.13) 


where  n^  is  the  number  density  of  the  electrons.  This  is  equivalent 
to  neglecting  the  shielding  by  the  ions  in  the  earlier  definition  of 
the  Debye  length  (2.3).  For  the  cases  of  attenuation  of  low  frequency 
oscillations  and  diffusion  across  a  magnetic  field  (i=s=l)  they  obtain 
a  screening  constant  between  (3.13)  and  (2.3)/  thus  indicating  that 
screening  by  ions  is  only  partially  effective.  For  the  cases  of  the  re- 
laocat  ion  of  an  anisotropic  distribution  of  ion  velocities,  and  of  thermal 
conductivity  across  a  strong  magnetic  field  (i=8=2)  they  again  obtain 
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tuw  same  result  as  in  (2.12)  with  the  effective  screening  constant 
between  (3*15)  and  (2.3)-  Thus,  in  these  cases  where  comparison  is 
possible,  ve  see  that  the  results  obtained  with  the  shielded  potential 
are  very  close  to  the  more  exact  results. 

A  difficulty  in  the  above  treatment  arises  in  verj"  dense  plasms 
when  the  Debye  length  d  turns  out  to  be  less  than  the  Interelectrcn 
distance  h  .  Thus  argon  at  1  atm  pressure  and  15000°K  has  only 
about  five  electrons  in  a  Debye  sphere.  In  this  case  the  concept  of  the 
Debye  length  as  a  shielding  distance  is  lost.  The  practice  adopted  in  CC and 
Cohen  et  al  [1950]  is  to  use  the  cut-off  procedure  with  the  distance  h 
rather  than  d  .  This  appears  to  be  the  only  alternative  with  the 
present  state  of  the  theory. 


3.2  Binary  Electron- Icn  Plasm 

This  case  has  already  been  considered  by  several  authors.  The  most 
accurate  and  detailed  work  seems  to  be  that  of  Landshoff  l19^9>  1951]  and 
Cohen  et  al  [1950]  and  Spitser  and  Harm  [1953]-  Landshoff  used  a  reduced 
form  of  the  Chapman-Snskog  method  out  to  the  fifth  approximation  (in 
the  notation  cf  this  report)  for  the  properties  in  the  absence  of  a 
magnetic  field.  Spitcer  and  co-workers  derived  their  res’ults  by  a  nu- 
metd.oal,  analysis  of  Eq.  (2.1)  with  a  Fokker- Planck  rather  mhan  a  Boltz¬ 
mann  collision  term.  Both  treatments  consider  a  plasm  with  no  mean 
velocity,  neutral  in  the  large,  and  with  no  pressure  gradients.  They 
assume  further  that  the  ions  have  a  Maucwellian  velocity  distribution 
euid  hence  do  not  contribute  to  the  heat  and  current  fluxes.  This  assump¬ 
tion  edso  reduces  the  number  of  terms  in  the  linear  equations  which 
result  from  the  Chapman-Enskog  mthod  used  ty  Landshoff. 


The  results  for  the  electron-singly  charged  ion  plasma  are  generally 
presented  as  the  coefficients  of  Eqs.  (2.J+7)  for  the  current  and  heat 

fluxes. 


J 


a? 


a 


ar 


-  pl  -  K’  ~ 

br 


(3. lit) 
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The  results  of  Spitzer  and  HSrm  for  these  coefficients  can  be  written 
as  follows: 


=  0.5908 


6 


and  by 


=  0.4155 


=  1.8904 


kfkT] 


Vm  e^ 


e  A 


(3.15) 


K'  =  2.2873 


Vm  e  S^n 
e 


As  before  the  true  thermal  conductivity  is  obtained  by  setting  the 
current  equal  to  zero.  They  obtained 


K  =  eK'  (=X) 


€  =  0.4189 


(3.16) 


Landshoff  [1951]  obtained  equations  of  the  same  form  as  the  above  with 
slightly  different  numerical  coefficients  in  his  fifth  approximation. 

His  results  will  be  presented  later. 

To  approximate  the  model  used  by  these  workers  with  the  computer 
program,  we  can  perform  the  calculations  for  various  ion/electron  mass 
ratios.  The  more  massive  the  ion,  the  smaller  will  be  its  average 
'velocity  and  the  less  will  be  its  contribution  to  the  transp^it  proper¬ 
ties.  Landshoff  neglected  the  ion  velocity  in  performing  the  integration 
over  the  impact  parameter.  Further,  to  agree  with  these  authors,  we  need 
to  consider  only  terms  of  dominant  order.  The  same  definition  of  A 
must  also  be  used,  since  this  has  been  taken  slightly  differently  in  the 
papers  referred  to  above. 

Examination  of  the  expressions  (2.14),  (2.15),  (2.26)  and  (2.34)  for 

Ip 

C.  .  ,  D.  ,  X'  and  q  in  conjunction  with  formulas  (2.42a),  (2.43a), 

IJ 


''In  a  binary  mixture  D  is  denoted  by  D.  following  MTGL. 

^  V  ^  J 
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(2.45a)  and  (2.46a)  shows  that  the  coeffit  ients  must  be  of  the  same 
form  as  those  of  Spitzer  and  Httrm  with  a  different  numerical  coefficient 
depending  on  the  level  of  approximation  and  on  the  ion/electron  mass 
ratio.  The  various  approximations  will  be  given  in  the  tables  to  follow 
after  division  by  some  standard  quantity.  Tables  1  and  2  give  the 
resu-lts  from  the  calculations  of  the  thermal  conductivities  K'  and  K 
after  division  by  the  first  approximation  to  the  electron  thermal  con¬ 
ductivity.  The  latter  may  be  determined  from 


w,  =  it 


>Jw  kT 


(3.17) 


2  (2  2) * 

with  wa  ^  as  given  in  (3-12).  The  results  of  the  other  authors. 

1/2 

who  essentially  neglect  terms  of  order  (m^/m^)  '  ,  are  given  in  the 

last  column  marked  »  .  The  level  of  the  approximation  in  the  notation 
of  this  report  is  given  in  the  first  column.  Since  the  results  of 
Spitzer  are  essentially  exact  to  dominant  order,  they  are  given  as  the 
infinite  approximation.  Ve  may  note  here  that  the  first  and  second 
approxxmations  of  Marshall  [1958]  agree  exactly  with  those  of  Landshoff 
[1951]  after  the  correction  by  Vaughn-Williams  and  Haas  [1961]  is  .made. 
The  first  approximation  as  derived  by  Chapman  [1954]  also  agrees  with 
these  results,  but  his  remark  that  the  higher  approximations  should  in¬ 
crease  his  resulcs  by  a  factor  of  about  1.4  is  not  in  agreement.  An 
explicit  formula  for  the  third  approximation  to  the  binary  plasma  thermal 
conductivity  has  also  been  derived  by  Inshennik  [I962].  For  v^^/a^=42.9 
he  says  that  his  formula  gives  [X] ^=0.741  [X^]  ,  which  also  agrees  with 
the  results  of  this  report. 


The  first  and  second  approxi-matiens  to  the  viscosity-,  after  division 
by  the  first  approximation  to  the  ion  viscosity,  are  given  ir.  Table  3- 
The  first  appx'Oximati>’n  to  the  ion  viscosity  is  given  by 


r  1  5  -v/TmikT 

*  ’’  6  2  (  2  2 )  * 


(3.18) 


2  (2  £ ) * 

with  the  average  cross-section  tt  a  given  in  equation  (3.12) 

lu  Tables  k  —  C  are  listed  the  other  coefficients,  p^te"  division  by  the 


corresponding  result  of  Spitzer  and  Harm. 
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TABLE  1.  HIGHER  APPROXmTIOI^'S  FOR  K*  FOR  THE  FULLY-IONIZED 

PLASMA  WITH  IONS  OF  SINGLE  CHARGE.  [K']  -^[X  ], 

nr  el 


II 

E 

IS36 

10^ 

10« 

10® 

m  =  1 

0.3242 

0.3127 

0.3042 

0.3033 

m  =  2 

1.4838 

1.4737 

1.4655 

1.4646 

m  =  3 

1.7461 

1.73^0 

1.7239 

1.7228 

E  =  4 

1.7591 

1.7471 

1.7389 

1.7358 

m  =-  5 

m  = 

X  Land she ff  [1951 

] 

*  Snit: 

;er  and  H*<rTi  [1953] 

0.3032X 

i.  J4-6ii-4x 
1.7227X 
I.7336X 
I.732IX 

1.7301* 


TAELE  2.  .HIGHER  APPROXIT-lATIOIiS  70  THE  THSRJ4AL  COrHDUCT: 

K  (=X)  FOR  A  FULLY- IONIZED  PLASMA  WITH  IONS 

SINGLE  CHARGE.  [K]  /[>.  ], 

mf  el 


1836 


0.3242 

c.32i;2 

0.7397 

0.7404 


:VITi‘ 

OF 


“  — 

m^/ 

TTl  ^ 

“e 

n 

=  1 

m 

II 

m 

=  3 

HI 

=  4 

m 

=  5 

a 

s  00 

0.3127 

0.312^ 

O.T269 

'.7275 


Key  as  for  Table  ; 


TABLE  3-  HIGHER  API 


10^ 

10^ 

0.3042 

0.3033 

0. 3032x 

C . 3042 

0.3033 

C.3032X 

0.7186 

0.7155 

0.7155X 

0.7172 

0.7161 

o.7i6c'x 

0.7183X 

0.7247* 

i 

HIGHER  AP.^ROXr-'J\TIONS  TO  THE  VlSCCSm'  FOR  A 
FULLY  IONIZED  PMSy.A.  V.TTH  IONS  OF  SINGLE  CHARGE. 


,  6 

0 

10 

10 

0.95&'' 

1.0851 


0.9981 

1.1485 


0.9998 

1.1514 


TABLE  k.  HICmR  APPROXIMATIONS  TO  TKE  HLECTRICAL  CONDUCTIVITY 
FOR  A  FULLY- IONIZED  PLASMA  WITH  IONS  OF  SINGLE 

CHARGE. 


1836 

4 

10 

CD 

m  »  1 

0.50^' 

Q,5o6h 

— 

0.5061+ 

0.5064 

0.5064x 

m  *  2 

0.9773 

0.9782 

0.978U 

0.9784 

0.9784x 

m  =  3 

0.9863 

0.9872 

0.987^ 

0.98/4 

. 987 4 X 

m  =»  4 

0.9933 

0.99^3 

0.99^5 

0.9945 

0.9934x 

m  =  5 

r . 7954x 

m  s  00 

1.0* 

[  Key  as  for  Table  1. 

TABLE  5.  HIGHER  APPROXLMATIONS  TO  ?  FOR  A  FJLLY- IONIZED 
PLASMA  WITH  IONS  OF  SINGLE  CHARGE.  [g]  / 


B9B 

1836 

lo** 

10*^ 

00 

0 

ao 

m  =  1 

0.3958 

0.3957 

0.3957 

0.3957 

0.3957X 

01  s  2 

i.0096 

1.01:2 

1.0103 

1.0103 

1.0103X 

m  =  3 

0.9903 

0.9909 

0.9910 

C . 9910 

0.99C^9x 

m  =s  4 

0.9983 

0.9989 

0.999C> 

0.9990 

0.9977X 

m  ®  5 

O.9987X 

n  =  00 

_ 1 

1.0* 

1  Key  as  for  Table  1. 

J 

_ 

TABLE  6.  HIGHEI.  APPROXLMATIONS  TO  a  FOR  A  FJLLY- IONIZED 
PLASMA  WITH  IONS  OF  SINGLE  CHARGE,  [a] 

m 


1837 

10^ 

10« 

108 

00 

m  =  1 

m  =  2 

m  =  3 

m  =  4 

m  =  5 

n  a  «o 

0 

1.1168 

0.9971 

1.0085 

0 

1.1181 

0.9980 

1 .0094 

0 

1.1134 

0.9982 

1.0096 

0 

1.1184 

0.9982 

1.0096 

0 

l.li84x 

0 . 9^1x 

i.oa/7x 

i.ro50x 

1.0* 

# 

Key  as  for  Table  1. 
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TABLE  7 .  higher  APPROXIMATIOfiS  TO  THE  THERMAL  CONDUCT'IVIIT 

X'  FOR  A  FULLY- lOraZED  PLASMA  VvTTH  IONS  OF 

SINGLE  CHARGE.  [X'  ]  /[X  ], 

itf  el 


m,/m  = 
i'  e 

1837 

10^ 

10« 

00 

m  =  1 

0.3242 

0.3127 

0 . 3042 

0.3033 

m  =  2 

0 . 6o64 

0.5891 

0.5868 

0.5859 

m  =  3 

0.9627 

0.9499 

0.9397 

0.9386 

m  =  4 

0.9668 

o.r‘  1 

0.9^37 

0.9^26 

In  all  of  the  results  we  note  that,  for  any  level  of  approximation, 
the  coefficients  are  fairly  independent  cf  the  ion/electron  mass  ratio. 

In  the  worst  case,  that  of  the  viscosity,  the  difference  is  only  about 

between  the  viscosity  of  a  hydrogen  plasma  and  the  hypothetical 

0 

plasma  whose  ion  has  10  times  the  i-  3S  of  the  electron.  The  corres¬ 
ponding  difference  for  the  thermal  conductivity  is  about  ,  and  less 
than  l']o  for  the  other  coefficients!  v;e  thus  have  rigorous  support  for 
the  simple  arg'mients  which  show  that  we  may  neglect  the  electrons  when 
ccmpi’.ting  viscosity,  and  neglect  the  ions  (but  not  the  electron-ion 
interaction)  when  computing  the  other  coefficients. 

We  can  new  examine  the  convergence  of  the  Chapman-Enskog  formulas 
fer  the  various  coefficients.  In  no  case  does  the  first  approximation 
give  a  satisfactory  answer.  The  second  approximation  to  a  and  g 
?.grees  within  a  few  per  cent  with  the  higher  approximations,  but  it  is 
necessary’  to  go  to  the  third  aprrexiaation  for  satisfactory  results  for 
n  and  a  .  Since  the  third  approximation  has  not  been  worked  out  for 
the  viscosity,  we  cannon  ray  with  certainty  how  much  the  second  approxi- 
matinn  differs  from  the  true  result.  From  study  of  the  rate  of  con¬ 
vergence  and  of  the  other  coefficients  ^t  would  apper.r  that  the  second 
approximaticr.  wculd  agree  within  a  few  per  cent  with  the  true  value. 

We  noted  earlier  that  the  f  rm  (2.23)  -'T  the  heat  flux  expression, 
which  was  used  in  this  section  to  afford  a  c , mparisen  with  the  results  of 
ether  authors,  is  not  the  most  convenient  form  for  computation  of  the 
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thermal  conductivity  of  multi-conponent  mixtures.  In  subsequent  calcu¬ 
lations  the  other  form  (2.23)  will  be  used,  sc  it  is  worthwhile  to 
compute  the  contribution  of  the  second  term  of  (2.24)  to  the  total 
thermal  conductivity  for  the  fully  ionized  gas.  This  is  done  in  Table 

where  [X']  is  given,  again  after  division  by  the  first  approximation 

m 

to  the  electron  thermal  conductivity.  V/e  note  that  the  thermal  conduc¬ 
tivity  would  be  seriously  overestimated  if  we  neglected  the  second  term 
of  (2.24).  This  term  is  proportional  to  the  thermal  diffusion  coeffi- 
ciects,  and  is  normally  taken  as  negligible.  They  are  still  quite 

small  in  this  case,  but  the  denominator  contains  the  term  m  la  ^  whicr. 

^  J 

is  also  small  when  both  1  and  ,J  refer  to  the  electron.  The  conclu¬ 
sion  is  that  the  contribution  of  the  second  term  of  (2.24),  negligible 
for  ordinary  gas  mixtures,  is  definitely  not  negligible  for  Ionized 
gases.  From  the  arguments  given  above,  we  can  probably  approximate  (2.2l) 
with  good  accuracy  for  the  pairtially  ionized  case  by 


X  *  X'  + 


OK 


n 


^  2 


(3.19) 


where  the  subscript  1  refers  to  the  electron.  In  the  special  case  of 
a  binary  ionized  mixture, 


=  D__m_ 
22  2 


but  ^22=0  so  and  only  E^2  ^21  non- vanishing.  But  also 

T  T 

D^=-D^  and  a^=n2  so 


X  =  X' 


E 

pk  12 


,Tx2 


(D^) 


n 


:Ei2.EjP  .  X' 


.  ££ 

n 


W2®12 


(3. 


20) 


This  approximation  for  the  case  of  more  than  2  species  will  be  examined 
in  a  subsequent  paper  for  partially  ionized  argon. 


-  41  - 


3 . 3  Mixtures  of  Isotopes;  Lorer.lzlan  Gas;  Quasi -Lorentz Ian  Gas. 


As  another  application  of  the  formulas  derived  in  Section  2,  we 
can  consider  three  special  gas  mixtures.  The  mixture  of  isotopes  of  the 
same  element,  which  will  supposedly  have  the  same  force  between  mole¬ 
cules  and  nearly  the  same  mass  is  the  first  case  to  be  considered.  Ve 
will  take  the  masseo  to  be  equal  so  the  results  will  be  the  same  as  for 
the  pure  gas,  except  .  .at  the  additional  coefficient  of  self  diffission 
will  be  considered  The  thermal  diffusion  coefficient  ■'/anishes  identi¬ 
cally  for  this  mixture.  For  the  Lorentzian  gas,  in  whlci.  one  compcnen 
of  the  binary  mixture  is  taken  tc  have  very  small  mass  and  concentra¬ 
tion,  there  are  exact  results  with  which  to  compare  the  Chapman- Enskog 
approximations.  Only  the  diffusion  coefficients  will  be  considered  here. 
The  last  mixtiua  to  be  considered  Will  be  that  of  the  quasi-Lorentzian 
gas,  a  binary  mixture  of  a  hea'/y  and  a  light  component,  the  former  of 
which  is  in  small  ccncentration.  For  this  mixture  Mason  [1957b]  dis¬ 
covered  the  exact  rsoult  for  the  binary  diffusion  coefficient.  He  has 
also  considered  these  mixtures  for  certain  of  the  inverse  power  poten¬ 
tials  usir.g  the  ChapmacA  and  Cc’^fling  formulation  for  the  binarj'’  mixture. 

For  the  inverse  power  potential 

cp  =  dr"^  (3.21) 

the  cross-sections  take  the  form  [S.3-*^] 
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Here-  is  a  pure  number  depending  only  on  6  .  The  Coulomb  case 

(Sal)  has  been  considered  earlier  in  a  somewhat  different  formulation. 
In  this  section  we  will  consider  only  terms  of  dominant  order  for  this 
case.  The  case  for  6=2  for  i=-l,2  has  been  worked  out  exactly  by 
ELiason,  Stogryn  and  Hirschfelder  [1956]  in  terms  of  sine  integrals.  We 
also  need  the  cross-sections  for  i=3>^  for  the  level  of  approximation 
considered  in  this  report.  These  can  be  worked  out  bv  the  same  method 
as  used  by  Eliason  et  al.  The  results  for  all  of  the  '(2)  are 


A^^^2)  -  -  f  [Si(27r)  -  2Si(7r)] 

A^^^(2)  =  -  J  [Si(l*7r)  -  2Si(27r)] 

(3.2M 

A^^^(2)  =  -  I  -  [Sl(2Tr)  -  2Si(7r)  ^  81(6^)  -  2Si(37r)l 

A^^^2)  =  -  ^  [SKStt)  -  4Si(27r)] 
where  Si(x)  is  the  sine  Integral 


Si(x) 


sin  Z 
Z 


dZ 


(3.25) 


These  expressions  may  be  evaluated  with  the  help  of  the  sine  integral 
tables  [F.W.A. ,19^0] .  The  results  a^e  given  in  Table  8. 


( i ) 

TABLES.  THE  QUAmTIES  k"  ’  {2)  OF  EQ.  (3-23)  FOR  i!=l,2,3,'+. 


1  =  1 

£  =  2 

II 

£  =  1+ 

II 

0.397601 

0.527  81+3 

0.712619 

0.812981 

The  case  for  6=3  has  also  been  considered  by  Eliason  et  al  with 

(i) 

an  approximate  technique  due  to  It'^tt- Smith.  The  A'  '^(3)  are  written 

ac  a  series 


.  1^3  . 


=  -2^1^  ^  Dp^  (nXsjrkr?)  <3-^) 

p=0  n=l 

with  the  Dp  as  pure  numbers  depending  only  on  5  .  They  have  been 
e'.'aluated  by  Eliason  et  al  for  p  up  to  4  .  Their  table,  which  evi¬ 
dently  contains  2Dp  rather  than  the  Dp  that  they  indicate,  is  re¬ 
produced  here  in  Table  9  after  dlvlsiou  by  2  .  The  values  determined 
for  A^^^(3)  are  given  in  Table  10  to  four  figures  though,  because  of 
the  slow  convergence  of  this  series,  they  are  probably  accurate  only 
cut  to  three  places. 


TABL^  9.  THE  COEFFICIS.NTS  Dp  FOR  THE  SERIES  (3.26). 


p  =  0 

P  =  1 

p  =  2 

p  =  3 

p  =  4 

Dp  = 

0.12719 

0.062435 

0.005095 

-C. 00141 

0.001585 

TABLE  10.  THE  QUANTITIES  A^*^3)  CF  EQ.  (3*23)  FOR  i=l,2,3,^. 


£  -  1 

i’  =  2 

i  =  3 

i.  =  4 

> 

It 

0.3116 

0.3535 

0.4724 

0.5034 

As  the  last  force  law  we  can  consider  the  gas  composed  of  hard- 
spheres  of  diameter  a  .  This  is  a  special  case  of  the  inverse  power 

2/  { )  2 

potential  where  we  take  the  limit  of  5-^  with  d  '  o  .  This 

case  may  be  treated  quite  simply  with  the  theory  of  Section  2  by  letting 
_]l  for  all  i  and  s  . 

The  results  for  the  mixture  of  isotopes  or  mechanically  simi.lar 
molecules  are  given  in  Tables  13.  and  12.  The  ratio  of  higher  to  the 
first  approximation  for  the  viscosity  and  thermal  conductivity  are  given 
in  Table  11  for  6=1,2, 3>4  and  »  .  The  case  of  5=4  is  known  as 
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Ilaxwelliaa  laolecules  and  gives  the  exact  result  in  the  first  approxim- 
tioa  to  all  coefficients.  This  case  was  not  computed  (though  it  would  he 
a  good  check  on  the  accuracy  of  th^  program)  but  is  merely  listed  for 
reference.  Because  of  the  exact  result  in  the  first  approximation  when 
5=4  we  might  expect  the  convergence  to  be  poorest  for  S=1  or  o=*  . 

We  see  from  the  table  that  5=1  do^s  indeed  give  the  slowest  convergence. 

TABLE  11.  HIGHER  APPROXIMATIONS  FOR  THE  THERT-IAL  CONDUCTIVITY  AND 
VISCOSITY  OF  A  MIXTURE  OF  ISOTOPES  FOR  THE  INVERSE 
POWER  POTENTIAL. 


[xyixl^  tnlj/tnljL 


m  =  2 

m  =  3 

1.02272 

1.02483 

1.01436 

4 

1.0 

1.0 

1.0 

3 

I.CO269 

1.00274 

1.00172 

2 

1.0250-0 

1.02508 

1.01579 

1 

1.25000 

1.26444 

1.15169 

,  1  II  ^  ■  .I..  1^ 

We  should  note  that,  because  of  the  vanishing  of  the  thermal  diffusion 
coefficient  for  this  case,  the  first  and  second  approximations  .must  be 
identical,  even  without  the  results  of  Muckenfuss  and  Curtiss  [1958] 
(See  Section.  2.3)-  For  this  reason  two  terms  in  the  Sonine  polynomial 
expansion  generally  are  called  the  first  approximation  to  the  theimal 
conductivity.  This  convention  is  followed  in  Table  11. 

In  Table  12  are  listed  the  self-diffusior  coefficients  for  the 
sane  values  of  8  as  well  as  for  others  which  have  been  worked  out  by 
Mason  [1957b] .  Here  again  the  slowest  convergence  is  found  for  5=1  . 


it  la  necessary  to  use  successively  smaller  ratios  of  these  quantities 
in  the  con5)utatlons  until  the  coefficients  appear  to  have  reached  their 
limiting  values.  Ratios  of  smaller  or 

larger  than  a  certain  size  will  cause  exponent  over-  or  underflow  in 
the  computer  when  certain  of  the  determinant  elements  are  evaluated, 
thus  leading  to  unknown  inaccuracies  in  thi"-  final  results.  Fortunately, 
convergence  of  the  formulas  occurred  for  ratios  which  the  program  could 
safely  handle. 

The  results  for  Lorentzlan  gas  are  presented  in  Tables  13  and  lU. 
Also  listed  are  the  mass  and  number  ratios  necessary  to  assure  that  the 
varloxis  approximations  had  converged  to  their  Lorertzian  values,  and  some 
results  taken  frrjm  y.ason  [1957b]. 


'TABLE  13.  HIGHER  APPROXIMATIONS  .’0  THE  DIFFUSION  COEFFICIENTS 
OF  A  LOREKTZIAN  GAS  FOR  INVERSE  POV/ER  POTENTIAL-3. 


■Mi 

■Mi 

m  =  3 

B 

exact 

— 

ay 

1.1068 

1.1165 

1.13177 

10-5 

10-5 

12 

1,048* 

1.05624 

0 

u 

10 

1.032* 

1.039* 

1.04528 

0 

0 

8 

1.023* 

1.027* 

1.03120 

0 

0 

6 

1.010* 

1.012*  ■ 

1.01373 

0 

0 

4 

1.0 

1.0 

1.0 

1.0 

0 

0 

3 

1.0119 

1.0131 

1.0135 

1.01373 

lo""^ 

1 

ON 

2 

1.1250 

1.1302 

1.1312 

1.13177 

lo""^ 

0 

1 

ON 

1 

. 

3.2500 

_ 

3.3906 

3.3945 

3.39531 

0 

1 

ON 

0 

1 

ON 

* 

Taken  from  Mason  [ 1957b]. 

TABLE  14.  HIGHER  APPROXBIATIONS  TO  THE  THEPJ^  DIFFUSION  RATIO 
lt)R  A  LORENTZIAN  GAS  FOR  IN\^RSE  POWER  POTENTIALS. 


L.  _ .. 

m  =  2 

m  =  3 

~  ■  1  — 

m  =  4 

V"a 

j 

00 

0.7692 

0.8941 

0.9388 

1C~^ 

10*^ 

j  12 

0.849‘» 

0.939* 

0.968* 

0 

0 

10 

0.865* 

0.947* 

0.972* 

0 

0 

!  8 

0.889* 

0.959* 

0.979* 

0 

‘0 

!  6 

0.928- 

0.975* 

0.998* 

0 

0 

i  4 

1.0 

1.0 

1.0 

0 

0 

i 

!  3 

\ 

i 

1.0588 

1.0135 

1.0054 

10  ‘ 

10"^ 

10-° 

1  2 

1.1111 

1.0133 

1.0043 

10  “ 

1  , 

1 

j. 

0.7692 

1.0138 

1.0012 

10-^ 

10-« 

*  Taken  from  Mason  [  1957^1  ]. 


We  note  that  the  error  in  the  first  approximation  to  the  diffusion 
coefficient  for  this  mixture  iS  very  large,  but  that  the  second  approxi¬ 
mation  is  quite  close  to  the  true  value.  The  rate  of  convergence  of 
the  thermal  diffusion  ratio  is  quite  slow  for  this  case.  Again  the 
inverse  single  power  potential  demonstrates  the  slowest  convergence. 
Fortunately  this  case  cannot  occur  in  practice  since  electrons  and  ions 
rust  be  present  in  nearly  equal  iensities  at  any  point  except  near  the 
tcuodaries  in  a  real  plasma.  Note  that  the  so-called  Lcrentzian  gas 
for  the  fully- ionized  plasma  [see  e.g.  Spitzer  and  Harm,  1953]  takes  the 
limit  but  retains  n^^n^  ■  Thus  the  ions  are  at  rest,  but 

the  electrons  and  ions  are  present  in  equal  densities. 

For  the  quasi-Lorentzian  mixture  ^^a3on  [195Tb]  has  discovered 
that  all  approximations  give  the  exact  result  for  the  binary  diffusion 
coefficient.  As  a  check  on  the  computer  program  we  verified  that  we 
could  duplicate  this  result  for  the  cases  6='», 3»2,1  .  Agreement  was 
achieved  out  to  better  than  six  significant  figures  with  the  mass  ratios 
listed  in  Table  15.  The  results  for  the  thermal  diffusion  ratio  are 


<il30  given  in  this  table.  No  exact 


result  exists  for  this  coefficient. 


TAHJ)  15. 


HIGHER  APPROXIMATIONS  TO 
FOR  THE  QUAoI-LORENTZIAN 

potential. 


THE  THERMAL  DIFFUSION  RATIO 
GAS  WITH  THE  INVERSE  P0LT:r 


m  *  3 

m  =  4 

V=>2 

00 

1.0455 

1.0515 

10"^ 

10^ 

12 

1.031* 

1.034* 

0 

0 

10 

1.028* 

1.031* 

0 

0 

8 

1.023* 

1.026* 

0 

0 

6 

1.016* 

1.017* 

0 

0 

k 

1.0 

1.0 

0 

0 

3 

0.9839 

0.9827 

10’^ 

loS 

2 

1.0909 

1.1098 

10"’^ 

10« 

1 

1.1360 

1.1742 

lO'"^ 

♦  Taken  from  Mason  [ 1957b] 

• 

We  note  again  the  slowest  convergence  for  the  Coulomb  potential,  though 
again  this  case  cannot  exist  in  practice.  The  rate  of  convergence  for 
this  mixture  is  considerably  better  t.han  for  the  Lorentoian  gas. 


k,  DISCUSSION  AND  CON'T.nSIONS 


In  this  report  the  usual  Chaprian-Enskog  formulation  of  the  transport 
properties  of  multicomponent  mixtures  has  been  extended  to  higher 
approximations.  The  expressions  derived  are,  strictly  speaking,  appli¬ 
cable  only  to  non-reacting  monatomic  gases  with  no  internal  degrees  of 
freedom.  However,  with  the  'xceptlcn  of  the  thermal  conductivity,  they 
can  probably  be  used  with  su'  •  .-ss  for  gases  in  which  reactions  are 
taking  place  and  even  for  polyatomic  gases.  The  corrections  necessary 
for  the  thermal  transport  in  monatomic  gases  were  given  in  Section  2.3. 

V7e  saw  in  Section  3-2  that  it  is  necessary  to  use  at  least  the 
third  approximation  to  obtain  accurate  results  for  the  thermal  conducti¬ 
vity  and  thermal  diffusion  coefficients  of  ion’’ oed  gases.  The  second 
approximations  to  the  viscosity  and  to  the  difliision  coefficient  of  the 
electrons  appear  to  give  quite  accurate  results.  The  theoretical  re¬ 
lations  are  quite  complicated  in  this  level  of  approximation  and  it 
would  be  worthwhile  to  derive  simpler  exrressicnw  which  approximate  the 
exact  results.  The  accuracy  of  any  approximate  relations  can  only  be 
Judged  by  the  agree.ment  ^or  lack  '’f  agreement)  with  the  accurate  results 
obtained  with  the  ex-nression?  given  here.  For  extensive  calculations  of 
the  properties  of  gases  it  would  appear  preferable  to  use  the  full  ex¬ 
pressions.  In  spite  of  complexity  of  the  formulas,  actual  computations 
’■roceed  quite  rapidly  on  a  typical  computer.  F^r  example,  the  coeffioie  "3  for 
argon  at  one  atmospaere  have  been  computed  from  to  22000*^K  in 

steps  of  1000°K  in  only  about  1-minute  run  time  on  the  IBM  7C9'"'w 
These  computations  were  performed  with  the  general  Balgol  program  men¬ 
tioned  earlier  and  the  time  in''luc]es  that  necessary  to  compute  seme  of 
the  avera.ge  cross-sections.  iGnised  argon  is  a  relatively  simple  mixture 
of  only  three  components,  and  somewhat  longer  time  wo\ild  be  required  for 
a  more  complicated  mixture. 

.Another  difficulty  in  usinp  the  general  expressions  arises  because 
not  ail  of  the  average  crcss-secticn  ’  necessary  for  the 

third  approximation  have  been  worked  out  for  the  potentials  of  inter¬ 
action  between  species  lo  more  complicated  mixtures.  Even  in  argon  we 
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■ire  United  to  the  third  approximation  because  the  integrals  necessary 
for  the  fourth  approximation  have  not  been  calculated  for  the  exponential 
IXJtentlal  which  is  used  for  the  atom-atom  and  atom- ion  elastic  inter¬ 
action.  Here  again  the  presence  of  high-speed  computers  should  make  it 
possible  to  carry  through  rapidly  the  cross-section  calculations  for 
various  poteatials.  This  is  not  done  for  the  argon  case  because  the 
third  approximation  is  felt  to  be  adequate.  In  cases  where  the  cross- 
sections  are  lacking  for  even  the  third  approximation  it  would  certainly 
be  worthwhile  to  set  up  a  general  program  to  compute  the  necessary 
cross-sections. 

One  additional  problem  not  completely  resolved  concerns  the  proper 

charged-particle-cross  section.  When  in  a  »  1  ,  then  the  Coulomb 

2 

potential  with  the  Debye  cut-off  is  adequate.  V/hen  a  »  1  but 
^  A  ^  1  then  a  more  accurate  method  must  be  used.  The  best  method  at 
present  appears  to  be  that  of  Kihara  and  Aono  [1963].  Because  of  the 
complexity  of  their  method,  the  simple  screened  potential  was  used  in 
this  report,  vi/here  comparison  was  possible,  the  results  agree  very  well 
with  those  of  Kihara,  Aono  and  Itikava  [1963]*  It  appears  that  no 
adequate  theory  has  yet  been  given  for  the  Coulomb  cross-section  when 
the  inter-electron  distance  becores  greater  than  the  Debye  length. 
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APPEIOIX 


la  this  appendix  relations  will  be  given  for  obtaining  the  bracket 
expressions  firoa  the  article  of  Mason  [I9?7a1*  Also  presented  are  tvj 
bracket  expressions  which  are  not  available  in  either  MITGL  or  Mason's 
article.  We  should  note  first  that  both  editions  of  OTGL  [1954,  1964] 
contain  a  misprint  in  Eq.  [7*A-5].  The  right  hand  side  of  this  relation 
should  be  multiplied  by  8  . 

Mason  [195’^a]  has  derived  the  higher  approximations  for  the  trans¬ 
port  properties  of  binary  gas  mixtures  with  the  formulation  of  Chapman 
acui  Cowling  [1952].  However,  he  does  not  use  exactly  the  same  expres¬ 
sions  that  they  present,  so  some  cross-checking  of  formulas  between  his 
article  and  Section  9*8  of  that  book  is  necessary  to  obtain  the  bracket 
expressions.  The  final  translation  formulas  are,  replacing  1  by  i 
and  2  by  J  in  the  expressions  of  Mason, 
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(m  >  0  >  m' 

(A.3) 


These  same  relations  may  also  be  used  to  derive  from  Mason  [1954]  some 
of  the  bracket  expressions  given  in  MTGL. 


T»/o  of  the  bracket  expressions  necessary  for  the  third  and  fourth 
approximations  in  the  MTGL  formulation  are  not  given  in  either  that  book 
or  in  Mason  [l957a]*  They  may  be  derived  from  relations  given  in  CC 
[p.  157,  Eq.  (1)];  see  also  the  remarks  in  Section  9.5  of  CC]  and  the 
bracket  expressions  already  obtained. 
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The  results  are 


!W  ®i 


2  35  (,(1,1)*  ^  o(l,2)*  + 


2  IJ 


*  6n<">3)' 


“i  __2  105  «(1,1)* 


(A. 4) 


-  56  - 


